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ABSTRACT 

Reconnection of the magnetic lines of force is considered in case the magnetic 
energy exceeds the rest energy of the matter. It is shown that the classical 
Sweet-Parker and Petschek models are generalized straightforwardly to this 
case and the reconnection rate may be estimated by substituting the Alfven 
velocity in the classical formulas by the speed of light. The outflow velocity 
in the Sweet-Parker configuration is mildly relativistic. In the Petschek con- 
figuration, the outflow velocity is ultrarclativistic whereas the angle between 
the slow shocks is very small. Due to the strong compression, the plasma 
outflow in the Petschek configuration may become strongly magnetized if the 
reconnecting fields are not exactly antiparallel. 
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1 INTRODUCTION 

In highly conducting plasmas, the magnetic energy is released by reconnection of the mag- 
netic lines of force near the singular lines where the magnetic field changes sign. This pro- 
cess has been intensively studied both in laboratory plasma devices and in space (see recent 
monographs by Biskamp (2000) and Priest & Forbes (2000)). Magnetic reconnection may 
occur in relativistic objects such as pulsars, magnetars, active galactic nuclei or gamma-ray 
bursts. In pulsars, the relativistic magnetic reconnection was proposed as a source of the 
high-energy emission (Lyubarskii 1996; Kirk, Skjaeraasen & Gallant 2002) and as the so- 
lution to the ex-problem (Coroniti 1990; Lyubarsky &; Kirk 2001; Kirk & Skjaeraasen 2003; 
Lyubarsky 2003). Similar models were also developed for the cosmological gamma-ray bursts 
(Drenkhahn 2002; Drenkhahn & Spruit 2002). The magnetic reconnection was evoked for 
explanation of the rapid variability observed in active galactic nuclei (Di Matteo 1998). The 



2 Y.E.Lyubarsky 

particle acceleration in the reconnection process was proposed to operate in radio jets (Ro- 
manova & Lovelace 1992; Birk, Crusius-Watzel & Lesch 2001 Jaroschek, Lesch & Treumann 
2004). The reconnection of the superstrong magnetic field is a key element of the widely rec- 
ognized model for the soft gamma-repeaters (Thompson & Duncan 1995; Lyutikov 2003). 
Therefore relativistic generalization of the classical reconnection models is in order. 

Blackman & Field (1994) considered kinematics of relativistic reconnection in the Sweet- 
Parker and Petschek configurations and concluded that due to the Lorentz contraction, the 
reconnection rate is significantly enhanced and may approach the speed of light. Lyutikov & 
Uzdensky (2002) confirmed this conclusion for the Sweet-Parker case. In both these works, 
the full energy and momentum balance was not considered, the authors imposed instead 
condition of incompressibility assuming that the proper density of the plasma remains con- 
stant. Particle acceleration in relativistic current sheets was studied both in the test particle 
approximation (Romanova & Lovelace 1992; Birk, Crusius-Watzel & Lesch 2001; Larrabee, 
Lovelace & Romanova 2003, Kirk 2004) and in two-dimensional PIC simulations (Zenitani 
& Hoshino 2001; Jaroschek et al 2004). 

Here we present generalization of the Sweet-Parker and Petschek models to the rel- 
ativistic case. In our analysis we will follow the classical approach just substituting the 
nonrelativistic expressions for the fluxes of the conserved quantities by the relativistic ones. 
It will be shown that the outflow velocity becomes relativistic if in the inflow region the 
magnetic energy density exceeds the plasma rest energy density. The classical formulas for 
the reconnection rate remain valid in this case if one substitutes the Alfven velocity by 
the speed of light. The relativistic Petschek reconnection differs qualitatively from the non- 
relativistic one in case the reconnecting fields are not exactly antiparallel. Due to very strong 
compression, a large field may be built up in the outflow region so that a significant fraction 
of the energy is ejected from the system as the Poynting flux. 

The article is organized as follows. The Sweet-Parker regime is considered in sect. 2. In 
sect. 3, the jump conditions at the relativistic slow shocks are obtained. These conditions 
are applied to the Petschek reconnection in sect. 4. Modification of the picture in case 
the reconnection fields are not strictly antiparallel are considered in sect. 5; the detailed 
consideration of this case is presented in Appendix. The results are summarized in sect. 6. 
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2 RELATIVISTIC SWEET-PARKER RECONNECTION 

In the Sweet-Parker configuration, the oppositely directed magnetic fields are separated by 
a current sheet with a small resistivity rj. Some anomalous resistivity is assumed taking 
into account that at high enough current density the sheet is known to be unstable to 
the growth of tearing mode and other instabilities. These instabilities are also developed 
in relativistic current sheets (Zelenyi & Krasnosel'skikh 1979; Lyutikov 2003). In the spirit 
of the classical reconnection models, we consider rj as a phenomenological parameter. The 
resistive dissipation of the magnetic energy within the sheet may be visualized as diffusion of 
the magnetic field towards the neutral plane where the oppositely directed fields annihilate. 
Outside the sheet, the magnetic field is frozen into the plasma and therefore magnetic 
diffusion brings the plasma into the sheet. The magnetic energy is converted into heat within 
the sheet and the thermal pressure thrusts the plasma out of the sheet through the sheet 
edges. The steady state is achieved when the plasma inflow is balanced by the outflow from 
the edges of the sheet. The reconnection rate may be roughly estimated from the integrated 
balance of the energy and momentum within the sheet. 

Let the current sheet be in the xz plane, the outer magnetic field be in the x direction 
(Fig.l). The sheet width is 2A, the sheet length is 21. The pressure equilibrium across the 
sheet implies 

where p is the plasma pressure within the sheet, B the magnetic field outside the sheet. 

The hot plasma in the sheet flows towards the edges of the sheet. The flow velocity in the 

sheet, v out , may be found from the momentum equation 

d 

fa (^OUt7oUt -P) = -JzBy, (2) 

where w is the enthalpy of the plasma. Here and thereafter the speed of light is taken to be 
unity. The current in the sheet is estimated from Ampere's law as 

Bo / q \ 

Substituting the x derivative by l/l, one gets 

^outTout ~P= ( 4 ) 
The vertical component of the magnetic field, B y , may be estimated from the flux conser- 
vation as By ~ ABq/1. Taking into account Eq.(l), one can see that the last term in the lhs 
of Eq.(4) is comparable with the rhs term. 
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Figure 1. Geometry of the Sweet-Parker reconnection. The current sheet is shaded. Magnetic field lines are shown by thin 
arrows, the plasma velocities by thick arrows. 



Let us consider the case when the magnetic energy density in the inflow region exceeds 
the plasma rest energy density. Then the plasma in the sheet is relativistically hot and one 
can take w = 4p. It follows immediately from Eq.(4) that the plasma may be accelerated 
to high Lorentz factors only if w significantly decreases towards the edges of the sheet. 
According to Eq.(l), p and, consequently, w are determined by the external magnetic field, 
which does not decrease significantly along the sheet. Therefore the plasma motion in the 
sheet is only mildly relativistic, 7 out ~ 1, u ou t ~ 1- Of course the relativistically hot plasma 
may be eventually accelerated to high Lorentz factors if the outer pressure considerably falls 
down. However the reconnection rate is determined by plasma motion in the region where 
the magnetic field is about its maximal value so that the pressure in this part of the sheet 
does not change much. 

Now let us consider the energy conservation. The energy influx is E z B /Att = v in B 2 /A7r, 
where t>i n is the inflow velocity determined just outside of the sheet where the magnetic field 
is frozen into the plasma. The energy outflow is wv out so that the energy balance is written 
as 

B 2 

v m -±l = Apv ont A. (5) 

47T 

Making use of Eq.(l), one obtains 

r, (6) 

^out f 

so that the flow is roughly incompressible (recall that all the above equations should be 
considered only as order of magnitude estimates). Therefore even if the reconnecting field 
are not exactly antiparallel and some small ^-component of the magnetic field is present in 
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the inflow region, the picture remains the same because B z remains small in the sheet and 
does not affect parameters of the outflowing plasma. 

In the steady state, V x E = 0, which implies dE z /dy = so that he electric field is the 
same within and outside the sheet. Within the sheet, E z obeys Ohm's law, which can be 
roughly written in the non- relativistic form, E z = r)j z , because the plasma motion within the 
sheet is only mildly relativistic. Outside the sheet, E z = v m B . Eliminating E z and making 
use of Eq.(3), one can estimate the inflow velocity as 

Eliminating A from Eqs. (6) and (7), one finds finally 

v m = S- 1 ' 2 - S = 4ttZ/t7 » 1. (8) 

So, contrary to what Blackmann & Field (1994) and Lyutikov & Uzdensky (2003) claimed, 
the inflow velocity in the relativistic Sweet-Parker regime remains much less than the speed 
of light. This is because the flow velocity in the sheet remains mildly relativistic so that the 
reconnection rate is not enhanced by the Lorentz contraction. 

3 JUMP CONDITIONS AT THE RELATIVISTIC SLOW SHOCK 

It was noticed by Petschek (1964) that the magnetic energy may be liberated not only in 
current sheets but also at slow shocks. Let us consider jump conditions at the slow shock in 
case the magnetic energy exceeds the plasma energy. Let the upstream flow be cold, W\ — p±, 
whereas the downstream flow relativistically hot, w 2 = ^p 2 - In the frame of reference where 
the shock is at rest and the upstream flow is perpendicular to the shock plane, conservation 

of the energy and momentum fluxes are written as 

2 . B a E t 2 B t2 E t 

PiliVi + = w 2 -f 2 v n2 + I (9) 

pi7iV + % = W ^ 2 A + H + P2-, (io) 

B n B t i 2 B n B t2 

—47- = - -£T- ( n ) 

Here the plasma density, p, and the enthalpy, w, are measured in the plasma rest frame 

whereas the electro-magnetic fields in the shock frame. Subscripts n and t refer to the 

normal and tangential components, correspondingly, and we take into account continuity of 

B n and E t . The last is written as 

E t = viB tl = v n2 B t2 - B n v t2 . (12) 
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Shocks arising in the Petchek model are close to the switch-off shocks for which B t2 = 0. 
Recently Komissarov (2003) demonstrated that, contrary to what was claimed before, such 
shocks are evolutionary not only in the non-relativistic but also in the relativistic case. For 
such shocks Eq.(12) yields 

v t 2 = (13) 

Substituting this relation into Eq.(ll) and eliminating w 2 with the aid of Eq.(9), one gets 

Bl{l-v l 

47^+^(1 

which means that the upstream flow moves with the Alfven velocity 1 as it should be. Intro 
ducing the magnetization parameter 



«? = ,,/T D2 ,r ..2v (14) 



one can see that 



(15) 



i>!=tan0 at a >> 1, (16) 

where 9 is the angle between the magnetic field and the shock plane. This assumes 9 < ir/4, 
which is fulfilled in the Petchek picture. In the opposite case the upstream velocity may be 
found from the full biquadratic equation (14); if B\ 3> 47rpi, 9 > it/ '4, the upstream flow is 
ultra-relativistic such that a is not large. Below we assume 9 < tt/A. 
Substituting Eq. (14) into Eq.(13) yields 



Vt2 = ~ C ° S 9 \ 1 I a 2n = " f 1 " o 1 2 J • ( 17 ) 

V 1 + a cos z 9 V 2a cos^ 9 J 
The last equality is obtained in the limit a ^> 1. Substituting Eqs. (16, 17) into Eqs. (9, 10), 

one can easily find in the same limit 

sin 9 

Vn2 = ~ 3^; 18 

2a cos 6 9 

72 = \[v cos 9; (19) 
B l cos2 

V, = (20) 

So the downstream flow is highly relativistic and directed at the angle ~ 1/a to the shock 
plane. 

Making use of the continuity equation, pi7it>i = P2l2V n 2, one finds the downstream 



1 In the plasma rest frame, the Alfven velocity is va = B' n j \J 4-n-p + B' 2 , where B' is the magnetic field in the plasma rest 
frame. One gets Eq.(14) taking into account that in the frame moving with the Alfven velocity, B n = B' n , B t = B[ (1 - ^)" 1/2 . 
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density 



p 2 =2/0i cos 2 9 



a 



(21) 



cos 29' 



So the downstream flow is highly compressed. The downstream temperature (or in fact the 
energy per particle) is easily found from the equation of state: 



where m is the average particle mass (m = m p /2 in the electron-proton plasma). Because 
T ^> m, the downstream flow may be loaded by pairs; in this case m becomes a function of 
the temperature. 

4 RELATIVISTIC PETSCHEK RECONNECTION 

In the Petschek (1964) picture, the current sheet is localized at — I < x < /, y = and 
pairs of slow shocks stem from the edges of the sheet as far as the outer boundary of the 
box at x — ±L (Fig.2). At the shocks, the x component of the magnetic field drops to 
zero and the magnetic tension pushes the plasma away from the reconnection region along 
the x axis. In response to rapid evacuation of the plasma from the reconnection region, the 
upstream plasma is sucked in toward the xz plane, together with the upstream field. The 
magnetic field in the inflow region is a small perturbation to a uniform horizontal field, 5 x. 
The plasma pressure and inertia are negligibly small at a >> 1; the electric force is also 
small because the upstream velocity is non-relativistic, as it can be checked a posteriori. 
Therefore the upstream field is potential, B = V$. According to the jump conditions, the 
plasma upstream of the slow shocks moves with the Alfven velocity, v- in = B y /B x m B y /B . 
It follows from the above analysis that the angle between the shocks is only about 1/a. 
Neglecting the inclination of the shocks, one should solve Laplace's equation in the upper 
half-plane with the boundary condition that the field becomes horizontal at large distances 
and that the normal component of the field be v m B§ between I and L at the x-axis and, by 
symmetry, — Vi n B between — L and — /. 

Following to Petschek (1964, see also Priest & Forbes 2000), one can write such a solution 

as 




(22) 





(23) 
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Figure 2. Pctschek rcconnection. The slow shocks are shown by thick lines; the rest elements are the same as in Fig.l 



Then the field at the origin is 



5,(0,0) = B 1 



2v in L 
ln- 

71 I 



(24) 



Taking into account that the mechanism chokes itself off when the field at the origin becomes 
too small, one can estimate a maximum reconnection rate by putting 5,(0,0) = 0.55 to 
give 

= 4 (25) 

The current sheet is described by the Sweet-Parker relation (8), which may be written as 
I = r] / (47ivf n ) . Substituting this into Eq.(25) and defining the Lundquist number via the 
external scale, S = 4nLc/ri, one gets finally 

7T 



™ *ms- (26) 

So the maximal reconnection rate may be estimated as about 0.1 of the speed of light. The 
reconnection rate does not approaches the speed of light, contrary to what was expected by 
Blackman & Field (1994), because in the outflow, the Lorentz contraction is compensated 
by a small angle between the slow shocks. According to Eq.(16), the inclination angle of the 



magnetic field in the inflow zone is 9 



0.1. 



5 THE RELATIVISTIC PETSCHEK RECONNECTION IN CASE B z + 

Downstream of the slow shocks, the plasma is highly compressed. Therefore if the reconnect- 
ing fields are not strictly antiparallel, the component of the magnetic field B z parallel to the 
current may become so large that the structure of the plasma outflow may be significantly 
affected. 
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Figure 3. Petschek reconnection in the oblique case. The current flows along the z axis; rotational discontinuities (RD) and 
slow shocks (SS) stem from the current sheet. A typical magnetic field line is shown by arrows. The trajectory of the point of 
intersection between the field line and the rotational discontinuity is shown by dotted line. 



Let us assume that some small B z \ = aB is presented in the inflow region (note that B z 
has the same sign in the upper and lower inflow regions). Then the upper and lower parts 
of the reconnected magnetic field line diverge in the z-direction. In this case, the shrinkage 
of the field line after the reconnection is accompanied by stretching in the ^-direction of 
the segment located within the outflow. Therefore the outflow should be separated from 
the inflow by rotational discontinuities where the magnetic field lines turn appropriately. It 
will be shown in Appendix that a pair of slow shocks should appear between the rotational 
discontinuities in order to satisfy the condition v z (y = 0) = 0. The structure of the flow is 
shown in Fig.3. In this section we only roughly estimate parameters of the outflow zone; the 
full solution is given in Appendix. 

Let ip r be the angle the rotational discontinuity makes with the xz plane. The magnetic 
field line shrinks with the speed of light therefore the point of intersection between the 
field line and the upper rotational discontinuity move in the z direction with the velocity a 
and in the y direction with the velocity ip r . At the bottom discontinuity, the corresponding 
velocities have opposite signs therefore the segment of the field line between the rotational 
discontinuities is inclined at the angle a/ip r to the xz plane. The y component of the magnetic 
field practically does not vary across the discontinuity. Therefore the ^-component of the 
field in the outflow may be estimated as 



It is taken into account in the last equality that in the inflow zone, the field line is inclined 



B 




(27) 
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at a small angle 9 to the xz plane (but 9 ^> ip r ). The estimate (27) may be from the 
conservation of the magnetic flux crossing a fluid contour, j>B ■ v x dl = 0. Choosing a 
rectangular contour in the xy plane with the rotational discontinuity as a diagonal, one 
finds v in B z i = B Zfint %p r . The inflow velocity is equal to the Alfven velocity (16) because 
the rotational discontinuity moves with respect to the plasma with the Alfven velocity (see 
Appendix). Then one immediately gets Eq.(27). 

The structure of the outflow is changed if the pressure of the B z field in the outflow zone 
becomes sufficient to compensate the outer pressure. In the rest frame of the outflow, the 
field is equal to B Z)OUt /^ out , where 7 out is the Lorentz factor of the plasma outflow. Taking 
into account that the upstream field is not Lorentz transformed because it is directed nearly 
parallel to the x axis, one concludes that the influence of B z is negligible if 

5 2 ,out < A)7out- (28) 
In this case the outflow parameters are described by Eqs.(18) and (28) like in the case 
B z = so that one can write 7 out = a/o 1 , ip = v n 2 = 9 /{2a). Taking this into account and 
eliminating B ZtOUt from Eqs.(27, 28), one can write the condition that B z in the outflow 
region is dynamically insignificant as 

a < (2^)-\ (29) 

In high-a plasmas, this condition is very restrictive; it may be violated at rather small a. 

In the opposite limit, the pressure and the energy in the outflow region are dominated by 
the magnetic field. The pressure equilibrium in the transverse direction yields (cf. Eq.(28)) 

B Z: out = -B 7out- (30) 
Increasing of the magnetic field at the rotational discontinuity does not contradict to the 
general principles because the field strength should remain the same only in the zero electric 
field frame. In the laboratory frame, the field grows significantly. The energy balance reads 
as 

VinB 2 = 5 2 2 out W. (31) 

Taking into account Eq.(27), one gets the half-angle between the rotational discontinuities 

as 

ijj r = 9a 2 . (32) 

and the Lorentz factor of the outflow as 

7out = l/a. (33) 
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One can see that at the condition opposite to that of Eq.(29), the opening angle of the 
outflow is larger whereas the outflow Lorentz factor is less than the corresponding quantities 
obtained at B z = 0. The transition to the B z = case occurs when a is about the rhs of 
Eq.(29). 

So even a small B z in the upstream flow significantly affects the structure of the outflow. 
At the condition reverse to that of Eq.(29) (but still a « 1), the strongly magnetized 
(-Sjz.out ~ Bo/a) plane jet of the angular width (32) is ejected with the Lorentz factor (33). 
In contrast with the B z = case, the energy flux is dominated by the Poynting flux. However 
the reconnection velocity is equal to the Alfven velocity in any case. Therefore the estimate 
(26) for the reconnection rate remains valid even if B z ^ 0. 



6 CONCLUSION 

The above analysis demonstrates that the classical reconnection models may be straight- 
forwardly generalized to the relativistic regime. This regime arises when the energy density 
in the inflow region is dominated by the magnetic field. Then the Alfven velocity is close 
to the speed of light. The reconnection rate in this case may be estimated by substituting 
v a = c in the classical formulas. The ejected plasma is relativistically hot so that the pair 
production is in principle possible. 

In the Sweet-Parker configuration, the plasma is pushed away with mildly relativistic 
velocities. In the Petschek configuration, the plasma is ejected with a large Lorentz factor 
within a very narrow angle. The qualitative difference with the non-relativistic reconnection 
arises if the reconnecting fields are not strictly antiparallel. At the non-relativistic slow 
shock, the compression ratio is finite so that if B z is small in the inflow region, it remains 
small and does not affect the outflow parameters. The plasma compression in the relativistic 
Petschek reconnection is very high therefore even at a rather small B lz the magnetization of 
the outflow may be very high so that most of the energy of the reconnecting fields is taken 
away by the Poynting flux. Therefore the relativistic Petschek reconnection should not be 
considered as a mechanism for the direct conversion of the magnetic energy into the plasma 
energy. However one can speculate that the ejected energy is eventually transferred to the 
plasma in the course of the flow expansion and/or development of MHD instabilities. 
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APPENDIX. RELATIVISTIC PETSCHEK RECONNECTION IN THE 
OBLIQUE CASE 

Let the magnetic field in the inflow zone have a nonzero component B z = aB . It was shown 
in Sect. 5, that one can neglect the obliquity at the condition (29). In this Appendix, the 
opposite case 

a 2 a > 1 (Al) 
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will be considered. The structure of the flow is shown in Fig. 3. The outflow is confined 
within two rotational discontinuities. Inside the outflow, two slow shocks occur. The rota- 
tional discontinuities and slow shocks are inclined to the xz plane by the angles ip r and ip s , 
correspondingly. At the rotational discontinuity, the magnetic field line turns whereas the 
plasma is pushed in the x direction. As we will see, the jump conditions require nonzero v z 
beyond the discontinuity whereas by symmetry, the plasma should move strictly parallel to 
the x axis at y = 0. Therefore a pair of slow shocks arises where the velocities are adjusted 
appropriately. 

Let us first find the jump conditions at the rotational discontinuities. Let us choose the 
coordinate system (£, n, () such that the upstream magnetic field lies in the £77 plane and the 
discontinuity in plane (Fig. 4). This coordinate system is turned with respect to that in 
Fig. 3 by the angle a around the y axis and by the angle tp r around the z axis. Taking into 
account that if) r <C a, 9, one can neglect the tangential component of the upstream velocity. 
Then the continuity of the tangential component of the electric field reads as 

E ( = VxBq = v v2 B i2 - v&B v \ (A2) 

E^ = = v (2 B v - v v2 B (2 . (A3) 

Here we take into account continuity of B v . 

It is known from the general theory that the proper density and entropy do not change 
at this sort of discontinuities. Employing this fact from the beginning, we need only three 
equations in order to find the rest of parameters. One can conveniently choose the continuity 
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of the flow, 

7i^i = 72%2, (-44) 
and the rjrj and r]( components of the momentum flux: 

Bl + B^-E^^Bl; (A5) 

Pllv^v^ ~ ^- (B V B (2 + E v2 E c ) = 0, (A6) 

47T 

where 

B B 

E v2 = B C2 vz 2 - B i2 v (2 = ~^ J1 v 1 . (A7) 

In the last equality, v^ 2 and v^ 2 were eliminated with the aid of Eqs.(A2) and (A3). 
Substituting Eqs.(A2) and (A7) into Eq.(A6) yields 

P72 2 ^ = ^(tan 2 0-u 1 2 ), (AS) 
where is the angle between the discontinuity and the upstream magnetic field. Then with 
account of Eq.(A4), one gets 

v\ = tan 2 e(l--), (A9) 



where a ^> 1 is defined by Eq.(15). So the upstream velocity is the same as in the aligned 
case therefore obliquity does not affect the reconnection rate. 

Now let us return to Eq.(A8) and eliminate, with the aid of Eqs.(A2) and (A3), v^ 2 and 
t>£ 2 from 7 2 . Then substituting v\ from Eq.(A9) yields the equation for v v2 : 

( B h + B h , l , i \ v 2 

{ Bl + + tan 2 9(1 - tan 2 9) J r > 2 

+2^v v2 - - = 0. (A10) 

In order to avoid cumbersome expressions, let us take into account that in the case of 
interest \B^ 2 \ \B^ 2 \ and, as follows from the preliminary estimates (see sect. 5) and 
will be confirmed below by rigorous evaluation, \B V /B^ 2 \ ~ 9a 9. Then the solution of 
Eq.(AlO) at the condition (Al) is written as 

iAU) 

The magnetic field downstream of the discontinuity is perpendicular to the x axis 2 (see Fig. 

2 Strictly speaking, the magnetic field line should be perpendicular to the x axis only at y = z = (by symmetry) that is only 
between the slow shocks. Above the shock, the magnetic field line lies in the plane set by the normal to the shock and the field 
line between the shock. Taking into account that tp s is much less than all other angles involved, one can safely assume that the 
field is perpendicular to the x axis already downstream of the rotational discontinuities. 
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3). Therefore 

% = -aB <2 . (A12) 

Then Eqs. (All), (A2) and (A3) yield at the condition (Al): 

B B 2 

v V 2 = v C2 = 2a; = -1 + 2a 2 . (413) 

Substituting Eqs.(A7), (A9) and (A12) into Eq.(A5),one finds 

= -aB (2 , (414) 

and then 

v v2 = -29a 2 ; 72 = (415) 

Now let us return to Fig. 3 and the coordinate system xyz. It follows form Eq.(Al4) 
that in the outflow, B z = Bo/a. Projecting v 2 onto the x axis, one finds v x2 = 1 — a 2 /2 so 
that the outflow moves in the x direction with the Lorentz factor 7 out = 1/a in accord with 
the estimate (33). By symmetry, the flow should move exactly in the x direction at y — 0. 
However projections of v 2 on the y and z axis are non zero. Therefore a pair of slow shocks 
arise in the outflow in order to adjust the velocities appropriately. 

Let us now find parameters of the slow shocks necessary to make the plasma move 
exactly in the x direction. Note that the quantities upstream of the shock are the same as 
downstream of the rotational discontinuity; we retain them with the index 2. The postshock 
quantities will be denoted by the index 3. The jump conditions may be conveniently found 
in the frame moving in the x direction with the Lorentz factor 7 out ; the physical quantities 
in this frame will be marked by prime (Fig. 5). In this frame, plasma moves only in the 
transverse direction with the Lorentz factor 

72 = 12/lout = l/2a, (416) 

whereas the shocks and rotational discontinuities move in the y direction with the non- 
relativistic velocities 

V's=i>slout = i>s/a and V/: = iJj r /a, (417) 

correspondingly. The velocity v 2 was calculated above with respect to the rotational discon- 
tinuity. With respect to the shock, the upstream normal velocity is 

<2 = < 2 + V; - V' a = -2a6 + (Vv - 1> a )/a. (A18) 
It follows from the estimates in sect. 5 and will be confirmed below that the velocities of the 
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Figure 5. Structure of the flow in the frame moving in the x direction with the Lorentz factor "font, view end on. The magnetic 
field line is shown by thin arrows, the rotational discontinuities (RD) and slow shocks (SS) by thick lines, velocities of the 
discontinuities by thick arrows. 

shock and rotational discontinuity are much less than the speed of light therefore the full 
Lorentz transformation is not necessary. 

Downstream of the shock, the plasma is at rest with respect to the "primed" frame of 
reference therefore with respect to the shock, tangential velocity is zero and the normal 
velocity is 

v' 3 = ~V: = -^s/oc. (A19) 



The angle between the upstream magnetic field and the shock is 
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(A20) 



Continuity of the tangential component of the electric field and of the normal component of 
the magnetic field take the form 

B' t2 (v' n2 - v' t2 <p 2 ) = B' t3 v' 3 ; (A21) 
B> t2 ip 2 = B' tm . {A22) 



Now the continuity of the energy and momentum fluxes may be written as 



P272H2 = 4 P37s^3 + 
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(A23) 
(A24) 
(A25) 



Inspection of Eq.(A25), with account of Eqs.(Al6, A18), shows that <£>3 — ip 2 ^2 at the 
condition (Al). Moreover one can safely neglect dynamical pressure as compared with the 
thermal and magnetic pressure in Eq.(A24). Then one gets 

(A26) 
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Substituting this relation into Eq.(A23) and eliminating P2l'<2 v 'n2 with the aid of Eq.(A25) 
(one can put v t2 = 1 and 73 = 1 there), one finds 

y? 2 = -Sv' 3 , (A27) 

which yields, with account of Eqs.(Al9, A20), 

ips = \ea 2 . (A28) 

Substituting Eq.(A27) into Eq.(A21), with account of Eq.(A20), yields 

v' n2 = -(4/3)^2 = -(4/3)0a. (29) 

Now one finds from Eqs.(A18), (A28) and (A29) that 

i> r = 9a 2 (A30) 

in accord with the estimate (32). 

One can find from the Lorentz transformation that B' n = B n /^ out = 9aB ; B' t2 = aB z2 = 
B . It follows from Eq.(A25) that 

1 - ? - A;- ^ 31 > 



= tt^tV- (32) 



Then Eq.(26) yields 

Bf 
I27ra 2 a 

Now one can find the ratio of the energy transferred away by the plasma and by the electro- 
magnetic field: 

B^/Att 3a V 1 ] 

So at the condition (Al), the energy of the reconnecting field is transferred away predomi- 
nantly in the form of the Poynting flux. 



